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ABSTRACT 


The  use  of  analytical  modeling  in  the  study  of  oceanic 
eddies  is  considered.  Limited  observational  data,  in 
combination  with  eddy  models,  can  be  used  to  obtain 
analytical  approximations  to  environmental  effects,  including 
current  and  temperature  perturbations,  throughout  the  eddy. 
Techniques  which  efficiently  use  discrete  measurements  are 
presented  to  accurately  specify  any  given  analytical  model 
containing  an  arbitrary  number  of  parameters  to  an  observed 
eddy.  Questions  of  unique  parameter  specification  and  data 
sufficiency  are  considered  for  various  data  types  and  amounts, 
using  a previously-derived  eddy  model.  Examples  with 
bathythermograph  data  are  presented,  in  which  eddy  size, 
strength,  and  center  position  are  to  be  determined.  AXBT 
data  is  emphasized,  and  an  investigation  is  made  of  the 
influence  of  the  number  of  such  instruments  on  the  accuracy 
of  parameter  estimates.  It  is  then  shown  how  data  obtained 
from  oceanographic  moorings  can  also  lead  to  specification  of 
eddy  drift  speed  and  direction.  In  both  the  bathythermograph 
and  mooring  examples,  it  is  demonstrated  that,  even  when 
the  type  of  data  available  leads  to  nonunique  parameter 
specification,  significant  information  can  be  obtained  about 
the  observed  eddy.  Results  in  this  paper  suggest  certain 
efficiencies  in  data  utilization  and  in  the  design  of 
subsequent  experiments. 
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1.  Introduction 

Oceanic  eddies  have  been  observed  to  be  responsible  for 
significant  environmental  perturbations  in  the  ocean's  static 
state.  Large  elevations  or  depressions  in  the  normally 
horizontal  isotherms  have  been  observed,  resulting  from  eddy 
rotational  currents,  and  the  translation  of  these  masses  of 
cold  or  warm  water  has  been  postulated  as  playing  a significant 
role  in  the  energy  balance  of  the  oceans  (Lai  and  Richardson, 

1977).  Investigations  of  individual  eddies  have  typically 
employed  extensive  deep  bathythermograph  readings,  to  depths 
of  several  thousand  meters  at  frequent  horizontal  locations  (see, 
for  example,  Fuglister,  1971) . Other  studies  have  also  incorporated 
current  data  obtained  from  moorings  (see  for  example,  Koshlyakov 
and  Grachev,  1974)  in  the  examination  of  eddies  in  mid-ocean 
regions.  In  such  experiments,  the  resulting  discrete  readings 
are  used  to  approximate  the  continuum  of  eddy  environmental 
effects  by  either  simple  smoothing  of  the  data,  or  the  more 
sophisticated  techniques  of  objective  analysis  proposed  by 
Bretherton,  Davis,  and  Fandry  (1976),  and  implemented  by  Freeland 
and  Gould  (1976).  However,  these  procedures  typically  require 
a high  density  of  readings  throughout  the  oceanic  region  containing 
the  eddy,  requiring  fairly  extensive  ship  time  and  equipment. 
Further,  results  from  these  procedures  are  restricted  to  regions 
within  the  eddy  where  data  was  acquired,  so  that  they  do  not  in 
general  determine  eddy  effects  at  all  horizontal  positions  and 
depths . 
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In  attempts  to  explain  and  describe  eddy  structure,  analytical 
research  has  produced  several  eddy  models.  Stern  (1975)  considered 
variational  principles  for  equilibrium  of  barotropic  eddies, 
with  possible  extensions  to  the  baroclinic  case.  Andrews  and 
Scully-Power  (1976)  proposed  a model  for  the  velocity  fields 
and  dynamic  heights  of  an  anticyclonic  eddy.  Flierl  discussed 
both  the  use  of  nonlinear  dynamics  to  study  eddies  linked  to 
the  barotropic  shear  flow  (Flierl,  1976) , and  linear  dynamics 
to  study  the  evolution  and  translation  of  eddy-like  disturbances 
(Flierl,  1977).  The  above  studies  did  not  attempt  to  compare 
theoretical  structures,  such  as  those  of  temperature  and  current/ 
to  those  of  observed  eddies.  However,  a recent  model  of  several 
of  the  authors  (Henrick,  Siegmann,  and  Jacobson,  1977)  was  shown 
to  give  satisfactory  agreement  to  current  and  temperature 
perturbations  observed  in  one  Gulf  Stream  ring.  The  model  is 
applicable  tc  a class  of  oceanic  eddies  satisfying  certain 
assmmptions.  We  remark  that  any  eddy  model  will  necessarily 
reflect  the  ultimate  objectives  of  the  modeler.  For  example, 
studies  of  biological  activity  in  Gulf  Stream  rings  (see,  for 
example,  Wiebe  et  al,  1976)  might  require  high  accuracy  in 
predicting  near-surface  effects,  while  long-range  acoustical 
propogation  studies  might  require  knowledge  of  eddy  influences 
at  great  depths.  We  note  also  that  no  one  model  could  be  expected 
to  appropriately  describe  all  eddies,  nor  could  it  completely 
describe  one  particular  eddy  from  its  nascent  to  terminal  stages. 

We  consider  here  the  problem  of  utilizing  analytical  models 
and  limited  data  in  the  description  of  properties  of  an  observed 
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eddy.  If  the  environmental  effects  of  an  eddy  could  be  approxi- 
mated to  the  accuracy  required  for  a particular  application,  then 
estimation  of  model  parameters,  using  the  model  and  limited  data, 
would  be  an  efficient  procedure.  Indeed,  such  an  approach 
might  require  significantly  less  extensive  measurements  than 
presently  used,  while  simultaneously  avoiding  the  known  hazards 
of  ocean  undersampling.  Of  course,  we  caution  against  the  use 
of  an  eddy  model  without  adequate  justification  or  verification, 
since  this  could  prejudice  the  ocean  measurements  and  lead  to 
grossly  inaccurate  descriptions.  In  this  paper,  we  also  make 
the  important  assumption  that  enough  evidence  is  available  to 
warrant  application  of  the  particular  model  chosen.  Our  purpose 
here  is  to  specify  the  types  and  amounts  of  data  necessary  for 
model-parameter  specification,  and  to  present  some  techniques 
to  implement  such  a procedure. 

In  Sec.  2,  a brief  review  of  a previously  derived  eddy 
model  is  presented.  The  assumptions  implicit  in  this  model  are 
outlined  and  the  model  parameters,  including  eddy  size,  depth, 
direction  and  speed  of  rotation,  and  drift  trajectory,  are 
presented.  The  problem  of  using  observational  data  to  specify 
the  parameters  implicit  in  any  eddy  model  is  discussed  in  Sec.  3. 
The  question  of  the  minimum  amount  of  required  data,  a necessary 
prerequisite  for  accurate  parameter  specification,  is  discussed. 
Then,  a technique  for  accurately  specifying  model  parameters, 
using  more  than  the  minimal  data, is  presented.  These  procedures 
are  next  applied  to  the  particular  model  of  Sec.  2.  Sec.  4 
discusses  the  use  of  temperature  and  current  data  in  selecting 
model  parameters,  and  the  processing  of  such  data  for  use  in 


the  parameter-specification  process.  Several  distinct  varieties 
of  such  data  are  then  used  to  specify  model  parameters,  for  the 
model  of  Sec.  2 as  an  example.  The  parameter-fitting  technique 
is  first  applied  in  Sec.  5 to  time-independent  data  obtained 
from  bathythermographs.  It  is  shown  that  if  data  through  the 
eddy  is  taken  in  a path  lying  in  a vertical  plane,  ambiguity 
in  the  position  of  the  eddy  center  can  result,  while  inclusion 
of  nonplanar  data  leads  to  unique  parameter  values.  Sec.  6 then 
investigates  the  use  of  oceanographic  moorings  to  study  both 
eddy  structure  and  translation.  Nonuniqueness  of  parameter 
values  results  if  only  a single  thermistor  string  is  used,  but  , 
inclusion  of  additional  thermistor  strings,  or  current  meters, 
results  in  unique  parameter  specification.  The  principal 
results  of  the  paper  are  summarized  in  Sec.  7. 

2.  An  eddy  model 

Using  the  basic  equations  of  fluid  motion,  Watson,  Siegmann, 
and  Jacobson  (1976)  developed  a method  for  analytically 
determining  approximate  environmental  effects  of  a class  of 
mesoscale  flows,  which  include  large  oceanic  eddies.  From 
these  results,  an  analytical  model  for  eddies  nearly  in  geostrophic 
and  hydrostatic  balance  was  derived  (Henrick,  Siegmann,  and 
Jacobson,  1977) , which  we  now  summarize  briefly  in  revised  form. 

Approximations  to  the  static  state  of  an  ocean,  with 
horizontal  surface  and  bottom  separated  by  distance  D,  can  be 
obtained  from  the  fluid  equations  by  setting  the  velocity  field 
equal  to  zero.  The  result  represents  a vertically-stratified 
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fluid  in  hydrostatic  balance,  with  static  density  p and 

s 

pressure  Pg  given  by 

2 

p (z)  = p_  + p_  g z/c  + TT(z)  (la) 

s o o o 

and 

z 

Pg(z)  = g / pg  (r)dr  . (lb) 

0 

In  (1) , z is  depth  measured  positively  from  the  ocean  surface, 

p is  the  surface  density,  g is  gravitational  acceleration,  c 
o o 

the  surface  sound  speed,  and  ir(z)  is  potential  density.  We 
select  it  to  have  the  observationally  reasonable  and  analytically 
convenient  form 

tt  (z)  = pQ  A[l-  (1+Bz) _1]  , (2) 

where  A and  B are  constants.  The  values  of  p Qr  cq,  A,  and  B 
are  chosen  by  fitting  oceanographic  data  from  the  area  under 
consideration.  For  example,  if  data  on  temperature  T is 
available,  a useful  state  equation  is  one  proposed  by  Eckart 
(1958) , relating  density  p and  pressure  P in  the  form 

F(p,P,S,T)  = 0 . (3) 

If  salinity  S is  assumed  constant  at  35  °/oo,  (3)  can  be  solved 
(Baer  and  Jacobson,  1974)  for  temperature  as  a function  of  p 
and  P.  Then,  values  of  pQ,  A,  and  B can  be  determined  by  fitting 
the  theoretical  temperature  structure  to  the  observed  data.  For 
instance,  by  utilizing  data  from  near  Eleuthera,  we  obtained 
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the  typical  values  of  pQ  = 1.02322  gm/cm3,  A = 0.00634,  and 
B = 2.312  km-1. 

To  determine  motion-induced  perturbations  to  the  static 
state,  a right-handed  Cartesian  coordinate  system  is  defined 
at  some  latitude  <f,  with  the  x-axis  oriented  positively  to  the 
east  and  y-axis  to  the  north.  Appropriate  scaling  of  the 
fluid  equations  leads  to  simplifications  valid  for  parameter 
ranges  that  may  include  mesoscale  eddies.  Motion-induced 
perturbations  p and  P to  the  static  pressure  and  density 
fields  are  sought,  so  that  the  total  pressure  and  density  may 
be  written  as 

P = Ps(z)  + P(x,y,z,t)  (4al 

and 

p = Ps(z)  + p(x,y,z,t),  (4bl 

where  t represents  time.  The  resulting  horizontal  current  V = (u,v) 
can  be  written  in  terms  of  a stream  function  as 


and 

v = u>  , 
rx 


(5a) 


(5b) 


while  density  and  pressure  perturbations  are  given  in  terms  of 

1 1>  by 
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p = (2ftp  g ^sin  <}> ) (5c) 

O Z 

and 

P = (2pQft  sin<f>)  \l>  , (5d) 

where  ft  represents  the  earth's  angular  velocity.  The  stream 
function  is  then  constrained  to  follow  a time-independent  quasi- 
geostrophic  potential  vorticity  equation,  which  is  equivalent 
to  the  conservation  of  potential  vorticity  following  a material 
particle.  The  quasi-geostrophic  balance  hypothesis  has  been 
used  by  other  eddy  modelers  (see,  for  example.  Stern,  1975), 
and  was  tested  in  the  MODE  experiment  (Bryden,  et  al,  1975), 
although  the  results  were  inconclusive. 

Solutions  that  exhibit  eddy-like  characteristics  are  then 
sought  to  the  potential  vorticity  equation.  Since  eddies  are 
often  observed  to  be  nearly  circular  features  drifting  through 
the  ocean  (see,  for  example,  Fuglister,  1971) , approximate 
solutions  are  considered  which  are  radially  symmetric  about  a 
central  axis  and  which  translate  with  time  along  some  path  with 
assumed  functional  form  [xQ(t),  yQ(t)].  Hence,  the  variable 

r2  = [x  - xQ(t)]2  + [y  - yQ(t)]2  (6a) 

is  introduced,  and  the  stream  function  is  written  as 

\p  (x,y,z,  t)  = x(r,z).  (6b) 

f 
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A reasonable  set  of  boundary  conditions  can  be  postulated 

from  observed  eddy  features.  For  example,  eddy  currents  and 

perturbation  pressure  and  density  are  required  to  vanish  at 

some  finite  distance  r from  the  eddy  center  and  at  some 

o 

depth  zQ  (<D) , below  which  eddy  effects  are  assumed  negligible. 

Other  assumptions , physically  motivated  from  surface  conditions 
of  unidirectional  rotational  currents,  finite  pressure,  and 
vanishing  vertical  velocity,  lead  to  additional  requirements 
for  the  solution.  Finally,  the  maximum  surface  current  speed 
is  normalized  to  the  velocity  scale  UQ.  An  approximate 
solution  to  the  boundary  value  problem  is  then  obtained  by 
superimposing  the  vertical  barotropic  and  first  baroclinic 
modal  solutions.  The  resulting  solution  (Henrick,  Siegmann, 
and  Jacobson,  1977)  is 

X (r , z)  = U0ro[Jo(alr/ro)_Jo(al)] [F(z)-F(zq)] , (7a). 

where 

Kx  = t Dg{2a;La2Cofi  sin<j)  [F  (0)  -F  (zQ)  ] }_1  , (7b) 

F (z ) = (l+Bz)"1/2{cos[(Y/2)ln(l+Bz)  ) 

+ (1/y) sin  [ (Y/2)ln(l+Bz)] } , (7c) 

and 

II 

Y = 2 ir/ln  (1+Bzq)  . (7d) 

| 
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In  (7)  , JQ  is  the  Bessel  function  of  the  first  kind  and  order 

zero,  = 3-83  is  the  first  zero  of  the  Bessel  function 

of  the  first  kind  and  order  one,  a ^ = .582  is  the  maximum  of 

J^,  and  In  is  the  natural  logarithm.  In  (7b)  the  plus (minus) 

sign  corresponds  to  an  anticyclonic  (cyclonic)  eddy  with 

clockwise  (counterclockwise)  rotational  velocity.  We  note 

that  both  the  separability  of  horizontal  and  vertical  dependence 

and  the  superposition  of  baroclinic  and  barotropic  modes, 

employed  in  the  development  of  this  model,  are  characteristic  of 

other  eddy  models  (see,  for  example,  Flierl,  1976). 

Resulting  velocity,  density,  and  pressure  fields  can  be 

determined  through  (5)-(7).  With  density  and  pressure  thus 

determined  at  any  point  in  the  eddy,  temperature  can  then  be 

obtained  from  the  inverted  Eckart  equation  (3)  when  S = 35  °/oo. 

In  order  to  determine  eddy-environmental  effects,  the  parameters 

characterizing  the  eddy  must  be  specified.  In  the  particular 

model  described  above,  we  need  to  choose  the  radius  rQ,  depth 

of  influence  z , maximum  current  speed  U , direction  of 
o o 

rotation  [the  sign  in  (7b)],  and  the  translational  trajectory 

described  by  x (t)  and  y (t) . 

o o 

We  note  that  the  above  model  has  been  shown  to  agree  fairly 
well  with  the  observed  features  of  one  Gulf  Stream  ring,  and 
that  it  will  be  used  as  an  example  in  subsequent  sections. 
However,  it  does  not  admit  such  possibly  significant  features  as 
radial  asymmetry,  salinity  variation,  or  near-surface  effects. 
Nonetheless,  these  limitations  do  not  affect  the  general 
objectives  of  this  paper,  since  any  eddy  model  which  contains 
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an  arbitrary  number  of  parameters  may  be  utilized  in  the 
procedures  to  be  discussed. 

3.  Parameter  determination 

Using  observational  data  to  specify  values  of  parameters 
in  an  appropriate  eddy  model  requires  the  matching  of  observed 
perturbations  in  the  static  ocean  state  with  eddy  effects 
predicted  by  the  model.  Wd^esmsider  here  the  problem  of  using 
eddy-induced  environmental  deviations  of  general  types  to 
accurately  select  model  parameters, 
a.  Minimum-data  parameter  specification 

The  matter  of  determining  how  little  data  of  a certain 
type  is  necessary  to  specify  a given  eddy  model  has  obvious 
theoretical  and  experimental  importance.  By  examining  the  minimum- 
data  problem,  much  can  be  learned  about  the  analytical  dependence 
of  the  model  on  its  parameters.  Moreover,  knoi  _dge  of  the 
minimum  data  set  is  a necessary  prerequisite  to  determining  the 
amount  of  data  sufficient  to  accurately  specify  model  parameters. 

Any  model  that  specifies  eddy-environmental  effects  depends 

on  position  and  time,  and  is  assumed  to  include  n parameters 

P = (p, »•••»?  )•  If  values  for  these  parameters  can  be  selected, 
l n 

eddy-environmental  effects  at  any  point  are  then  known  to  the 
accuracy  of  the  model.  For  fixed  location  and  time,  the  eddy 
model  can  be  regarded  as  specifying  eddy  effects  as  nonlinear 
functions  of  the  n parameters.  Suppose  a deviation  A from  the 
ocean's  static  state  is  observed  at  a particular  position  and 
time  (x^,y^,z^,t^) . A corresponding  relation  of  the  form 


for  the  deviation  M predicted  by  the  model  equations  can  be 
obtained,  where  the  components  p^  of  are  viewed  as  independent 
variables.  It  is  tempting  to  suppose  that  given  a sufficient 
number  (>n)  of  observations,  a unique  choice  of  the  n para- 
meters could  be  determined  by  solving  equations  (8)  with  each 
M replaced  by  the  corresponding  A at  n distinct  positions  or 
times.  However,  due  to  the  functional  forms  of  (8) , it  is 
possible  that  these  measurements  may  not  be  taken  at  arbitrary 
points  in  space  or  time.  Moreover,  unique  parameter  values 
may  not  result,  even  if  a proper  sampling  is  performed.  Thus, 
the  problem  of  parameter  determination  requires  further 
consideration. 

We  suppose  here  that  the  model  equations  at  any  (x,y,z,t) 
are  separable  into  a product  of  vertical  and  horizontal  functions, 
each  with  distinct  parameters: 

M = fR(x,y,t;  jpjj)  fy(z  t;  _pv)  . (9) 

In  (9)  , p„  represents  the  vector  of  n„  parameters  p^^  , 

'M~.il  n n 

1 < i < n„,  associated  with  the  description  of  the  horizontal 
structure  and  translation  of  the  eddy.  Similarly,  pv  is  the 
vector  of  nv  parameters  p^^ , 1 £ i £ ny,  which  describe  the 
eddy  vertical  structure.  Such  separability  is  characteristic 
of  many  eddy  models,  and  appears  justifiable  from  observations 
of  near-uniform  decay  of  eddies  with  depth.  However,  if  such 
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separability  were  not  permitted  in  the  model,  it  would  be 
necessary  to  choose  all  the  parameters  at  once. 

Suppose  that  deviations  in  the  ocean's  static  state 
are  observed  at  a particular  fixed  horizontal  position 
and  time  (x^y^t^)  and  at  ny  + 1 different  depths  z ^ , 

1 < j riy  + 1.  The  vertical  parameters  can  be  determined 
by  equating  observed  and  predicted  environmental  deviations 
for  each  j , giving 

= fH(xi,yi,ti'PH)  fv(zij)  »ti;^v)  ' 1-^-nV+1  * tl0) 

The  horizontal  function  f„,  which  is  constant  at  all  depths, 

n 

can  then  be  eliminated  from  these  equations  by  solving  for 
fH  at  a single  depth,  and  substituting  into  the  remaining  n^ 
equations.  This  results  in  n^  equations  of  the  form 

A^Vz*1*,^;  pv)  - A^  fv(z^  ,ti;  J?v)  = 0,  2<j<ryH  , 

(ID 

where,  without  loss  of  generality,  we  have  chosen  the  first 
equation  to  eliminate  horizontal  dependence.  Approximate 
solutions  to  these  equations  for  the  n^  vertical  parameters 
can  then  be  obtained  numerically,  and  the  vertical  structure 
of  the  eddy  completely  specified. 

Knowledge  of  the  vertical  structure  can  be  used  to 
eliminate  vertical  dependence  of  measurements  at  any  point  in 
the  eddy.  If  a measurement  A^^  is  taken  at  any  horizontal 
position  (x^»y^)  and  time  t^,  and  depth  zj^ , depth  dependence 
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can  be  eliminated  by  defining  the  vertically-scaled  environ- 
mental deviation 

= A^^/fy  (z^  »t^;  Py)  » 1 £ j £m  , (12) 

where  typically  measurements  would  be  available  at  m depths 
^ , 1 £ j £ m.  The  magnitude  of  eddy  effects  at  (x-^/y^ft^) 
can  then  be  determined  by 

m 

l1  = (1/m)  ^ A^j)  , (13) 

j=l 

where  A^  is  the  average  vertically-scaled  environmental 
deviation.  This  averaging  procedure  avoids  propagation  of 
errors  in  data  at  any  one  depth,  and  hence  allows  more  accurate 
measurement  of  the  strength  of  eddy  effects  at  (x^»y^ft^). 

We  note  that  if  the  vertical  model  structure  fits  the  observed 
vertical  structure  exactly,  then  the  terms  in  sum  (13)  are 
identical. 

The  horizontal  parameters  can  then  be  chosen  using  the 
vertically  independent  measurements  from  (13) . Given  observations 
at  n„  distinct  horizontal  positions  and/or  times  (x. ,y. , t. ), 

n ill 

1 < i < n„,  where  at  least  one  of  the  x. ,y. , or  t-  are  distinct 
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for  each  i,  the  average  vertically-scaled  environmental  deviation 
A^  can  be  determined  at  each  position  as  in  (13) . Each  A^  can 
then  be  equated  with  the  vertically^scaled  deviation  predicted 
by  the  model,  resulting  in  nfl  equations  in  the  form 
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(14) 


for  the  n„  unknown  parameters.  Approximate  solutions  for  the 
horizontal  parameters  can  then  be  obtained  numerically, 
resulting  in  specification  of  horizontal  eddy  structure. 

However,  here  and  in  the  vertical  problem  (11)  it  is  possible 
that  exact  solutions  might  not  exist,  since  any  given  model 
may  not  exactly  match  observed  effects,  which  themselves  contain 
experimental  inaccuracies.  This  problem  can  be  alleviated  by 
using  procedures  outlined  in  the  following  subsection.  It  is 
also  possible  that,  as  a result  of  the  functional  forms  of  the 
model  equations,  nonunique  parameter  values  may  exist  for  certain 
types  of  data.  Examples  of  this  situation  will  be  given  in 
Secs.  5 and  6.  The  minimum  amount  of  data  required  to  completely 
specify  the  eddy  model  is  riy+1  measurements  at  different  depths 
(at  the  same  horizontal  position  and  time)  and  n„  measurements 

xl 

at  different  horizontal  locations  and/or  times. 

b.  Accurate  parameter  specification 

To  obtain  a more  accurate  analytical  approximation  to  an 
observed  eddy,  it  is  very  desirable  to  utilize  more  than  the 
minimum  required  data  set,  since  no  model  can  be  expected  to 
give  an  exact  fit  to  the  observed  process.  In  addition,  some 
of  the  inaccuracies  in  ocean  measurements  can  be  decreased  by 
employing  additional  measurements. 

The  problem  of  generating  more  accurate  choices  of  model 
parameters  can  be  viewed  as  overspecifying  systems  (11)  and 
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(14)  by  using  more  points  (and  hence  more  equations)  than 
parameters.  As  discussed  previously,  it  is  appropriate  to 
first  choose  the  vertical  parameters,  and  then  the  nfl 
horizontal  parameters.  If  horizontal  position  and  time  are 
fixed  at  (x^/y^t^)  and  observations  A^-^  are  obtained  from 
mv(>nv+l)  depths  z^^ , the  vertical  parameters  can  be  varied 
to  obtain  a reasonable  fit  at  all  depths.  The  choice  of  these 
parameters  could  be  made,  for  example,  by  minimizing  the  sum 
of  the  squares  of  error  of  the  fit.  However,  we  have  found 
that,  to  avoid  biasing  the  fit  to  the  larger  near-surface 
observations,  a more  accurate  agreement  with  observed  eddy 
structure  is  obtained  by  minimizing  the  sum  of  the  relative 
squares  of  the  error: 

"V  .2.1 

Min  £ { [Mfv(z^),t1;  Pv)  - A^h  /( A^*)2}  , (15) 

j-1 

where  M = fy  (x^ ' y^ ' t^ .We  note  that  the  magnitude  of  the 
horizontal  structure  M is  varied  since  eliminating  this  quantity, 
using  one  of  the  measurements  as  in  (11) , can  lead  to  significant 
errors  if  an  inaccurate  measurement  is  used  for  the  elimination. 
Problem  (15)  can  then  be  solved,  for  example  by  numerical 
techniques  such  as  a modified  Levenburg-Marquardt  algorithm 
(Brown  and  Dennis,  1972),  or  by  statistical  nonlinear  regression 
techniques,  to  minimize  the  error  in  interpolating  between  data 
points.  This  results  in  specifying  the  vertical  structure.  We 
note  that  for  accurate  specification  of  eddy  vertical  structure, 
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we  have  found  from  our  numerical  experience  that  it  is  necessary 

to  utilize  readings  to  a depth  of  750  m or  more. 

The  choice  of  the  n horizontal  parameters  requires 

H 

mH  ^ > nH^  measurements  at  distinct  horizontal  positions  and/or 
times.  At  each  position  and  time  (x^,y^,t^) , the  average 
vertically-scaled  deviation  is  computed  as  in  (13).  The 
resulting  minimization  problem  to  be  solved  in  order  to  find 
the  horizontal  parameters  is 


Min 

PH 


£ {l5i  ' £H<*i'yi'VjH,l2/<V2  } ' 

i=l 


(16) 


where  relative  least  squares  is  again  advantageous  to  avoid 
biasing  larger  terms.  We  note  that  more  accurate  parameter 
fits  will  be  obtained  by  using  a wide  distribution  of  points 
throughout  the  eddy,  rather  than  a narrow  grouping  of  points 
in  one  area,  such  as  near  the  eddy  edge.  Further,  because  of 
the  ordinarily  large  length  scales  of  eddies,  the  points  should 
be  separated  by  at  least  10  km  to  insure  distinct  readings; 
equivalently,  if  data  is  taken  at  the  same  horizontal  location, 
but  at  different  times,  intervals  of  several  days  between  points 
used  in  (16)  is  sufficient,  assuming  typical  eddy  drift  speeds 
of  3-5  km/day. 

Thus,  knowledge  of  the  minimum  required  data  facilitates 
overspecification  of  the  problems  for  the  vertical  and  horizontal 
parameters.  Such  overspecification  reduces  the  effects  of  errors 
in  measurements  and  deficiencies  in  the  model  to  give  a more 
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accurate  analytical  approximation  to  the  eddy  structure.  If 
the  environmental  effects  of  an  eddy  are  first  observed  at 
mv(>nv+l)  depths  at  the  same  horizontal  position  and  time,  the 
riy  vertical  parameters,  and  hence  the  eddy  vertical  structure, 
may  be  specified  by  solving  (15).  Then, given  m^  (>nH)  measure- 
ments of  environmental  deviations  at  distinct  horizontal 
positions  and/or  times,  the  remaining  n^  horizontal  parameters 
can  be  determined  by  solving  (16).  Thus, the  eddy  structure  is 
specified  at  all  locations  and,  if  the  original  data  is  taken 
at  distinct  times,  at  all  times. 

c.  Specialization  to  one  eddy  model 

We  conclude  this  section  by  specializing  the  discussion  of 
previous  subsections  to  the  authors'  model  of  Sec.  2.  The 
parameters  of  this  eddy  model  are  the  direction  of  rotation 
(specified  by  the  sign  in  (7b)),  depth  of  influence  zq,  radius 
rQ,  maximum  rotational  current  speed  UQ,  and  nD  parameters  for 
the  form  assumed  for  the  eddy  trajectory.  For  example,  if  a 
constant  drift  velocity  is  assumed,  the  position  of  the  eddy 
center  at  time  t is  given  by 

*xo(t)'  yo(t)]  =[Xo  + V'  Yo  + VDt]  ' (17) 

where  (Xq,Yo)  is  position  of  the  eddy  center  at  time  t = 0, 

and  = (UD'VD)  is  tJie  horizontal  drift  velocity.  Thus,  there 
are  nQ  = 4 drift  parameters  in  this  case.  Simpler  or  more  general 
forms  for  the  drift  trajectory  may  be  assumed,  but  nD  2 since 

II 
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two  parameters  representing  the  initial  position  of  the  eddy 
center  will  always  be  present.  We  observe  that  the  direction 
of  rotation  can  be  chosen  unambiguously  from  the  sign  of  the 
perturbation  temperature,  so  that  there  are  nH  = nD+2  horizontal 
parameters.  Moreover,  n^  = 1 since  the  only  vertical  parameter 
is  zQ.  For  any  parameter  values,  deviations  in  pressure  and 
density  structure,  and  current  velocities,  are  obtained  from 
(5),  (6),  and  (7)  as 

P = (2 p fisin  <J> ) U r K,  [J  (a.r/r  ) - J (a1)][F(z)  - F (z  ) ] , (18a) 

O OOIOXO  Ol  o 

p = (2pong-1sin  ^^r^  [ JQ  (c^r/r^  - J^c^)  )F'  (z)  , (18b) 

and 

(u , v ) = (air/ro)  IF  (z)  -F  (zQ)  ] [-y+yQ  (t)  ,x-xQ  (t)  J/r,  (18c) 

The  resulting  temperature  is  then  obtained  from  (l)-(3).  Thus, 
this  model  is  separable  in  the  sense  of  subsection  a. 

The  sole  vertical  parameter  zQ  may  first  be  determined  by 
measuring  environmental  deviations  at  a fixed  horizontal  position 
and  time  and  at  (>2)  different  depths,  solving  problem  (15) 
with  Pv^  = zq,  and  selecting  the  appropriate  vertical  function 
from  (8).  The  horizontal  parameters  can  then  be  obtained  using 
mH  ^>nH^  readings  at  distinct  horizontal  locations  or  times, 
and  then  solving  (16)  with  the  appropriate  parameters  and 
horizontal  structure  from  (18).  We  shall  illustrate  this  pro- 
cedure in  following  sections. 
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4.  Typical  observational  data 

We  consider  here  the  use  of  temperature  and  current 
measurements  to  select  model  parameters.  Our  purpose  is  to 
consider  the  efficient  use  of  such  data,  in  the  sense  that 
problems  (15)  and  (16)  are  numerically  stable  and  that  all  the 
information  implicit  in  the  data  is  utilized. 

a.  Temperature  measurements 

The  most  commonly-observed  environmental  effect  of  oceanic 
eddies  is  the  large  temperature  perturbations  induced  by  their 
strong  rotational  currents.  Each  eddy  produces  a characteristic 
temperature  perturbation,  a fact  which  has  been  used  in  the 
tracking  of  eddies  over  extended  periods  of  time  (see,  for 
example,  Parker,  1971) . Consequently,  it  is  of  obvious  importance 
to  consider  specification  of  model  parameters  using  temperature 
observations.  In  the  process  of  constructing  an  eddy  model, 
however,  temperature  is  typically  not  a primary  quantity. 

Solutions  to  the  fluid  equations  are  usually  in  terms  of  pressure, 
density,  and  current.  Resulting  temperature  effects  can  then 
be  determined  by  using  a state  equation  such  as  (3) , or  tabulated 
values  of  sea  water  properties.  Unfortunately,  this  determination 
is  typically  numerically  unstable,  in  that  small  changes  of 
perturbation  density  lead  to  large  temperature  variations,  making 
accurate  selection  of  parameters  difficult.  To  avoid  this  problem, 
we  propose  preprocessing  temperature  measurements  so  that  they 
may  be  used  effectively  in  parameter  selection. 


Eddy  temperature  effects  result  from  perturbations  in  the 
static  density,  salinity,  and  pressure  profiles.  By  assuming 
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small  perturbations  in  these  quantities  from  static  values  as 
in  (1),  and  expanding  (3)  in  the  perturbed  quantities,  it  can 
be  shown  that  the  resulting  temperature  perturbations  are 
primarily  functions  of  the  perturbation  density.  Physically, 
this  results  from  the  fact  that  eddy-induced  salinity  variations 
are  relatively  small  (Fuglister,  1972) , and  that  pressure 
perturbations  are  small  in  comparison  to  the  large  hydro- 
static pressure.  We  may  then  approximate  the  density  by 

p = E[T,Ps(z)]  , ' (19) 

where  (19)  represents  either  a form  of  sx ate  equation  (3)  with 

constant  salinity  or  else  a tabulated  density  value,  and  P (z) 

s 

is  the  static  pressure  in  (lb).  We  note  that  E(T,P)  is  a well- 
conditioned  function  of  temperature,  in  that  small  errors  in 
temperature  will  cause  smaller  errors  in  density.  Moreover,  the 
induced  error  in  approximating  density  in  (19)  can  be  shown  to 
be  less  than  one  percent  of  the  eddy-induced  perturbation  density. 

The  static  state  of  an  ocean  area  under  consideration  is 
typically  known  from  long-term  observations  or  archival  sources. 
Eddy  density  perturbations  Ap  may  then  be  approximated  from  the 
observed  temperature  perturbations  AT  by  the  relation 

Ap  = E (T  + AT,  P ) - E (T  , P ) , (20) 

o S S 5 

where  Tg  is  either  the  known  static  temperature  or  the  theoretical 
static  temperature  derived  from  (1)  - (3) . Equation  (20)  gives 
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perturbation  density  sufficiently  accurate  for  selection  of 
model  parameters.  Theoretical  density  predictions  from  the 
eddy  model  can  then  be  used  in  comparison  with  those  of  (20) 
in  the  procedure  of  Sec.  3. 

b.  Current  measurements 

In  recent  years,  current  measurements  have  become  increasingly 
important  in  studying  large-scale  oceanographic  phenomena  (see, 
for  example,  Treshinikov  et  al,  1977) . Clearly,  significant 
information  about  eddies  can  be  acquired  through  the  use  of 
i current  meters  (Koslyakov  and  Grachev,  1973). 

A single  current  meter  provides  both  the  magnitude  |v|  and 
the  direction  0 of  the  observed  current,  where  we  take  the 
i angle  0 to  be  measured  positively  counterclockwise  from  an 

eastward  latitudinal  parallel.  As  with  density  measurements, 

\ observations  of  current  magnitude  can  be  used  with  predictions 

of  current  speeds  from  the  model  equations  directly  in  the  pro- 
* cedure  of  Sec.  3.  On  the  other  hand,  current-direction  readings 

| provide  additional  information  about  the  position  of  the  eddy 

center,  i.e.,  xQ(t)  and  yQ(t).  How  this  information  is  exploited 
j depends  on  the  eddy  model.  In  one  with  assumed  circular  stream- 

lines, such  as  the  model  of  Sec.  2,  the  unit  tangent  vector 
'!  (cos  0,  sin  0)  to  a streamline  at  a given  point  will  be 

perpendicular  to  a line  through  the  eddy  center.  Thus,  if  a 
current  direction  0^  is  observed  at  the  horizontal  position 

(x. ,y . ) at  time  t . , the  eddy  center  is  known  to  lie  on  the  line 
11  1 

[i 

n 
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[yQ(ti)-yi]  = [xq  (ti) -x^^]  tan  (ir/2-6^^)  . (21) 

Equation  (21)  is  a relation  between  the  coordinates  (xq  (t^) , yQ (t^) ) 
of  the  eddy  center  at  the  time  of  the  measurement.  If  current 
direction  is  measured  at  two  or  more  points  simultaneously,  and 
if  the  points  do  not  lie  on  the  same  eddy  diameter,  then  xQ(t^) 
and  y (t±)  are  uniquely  determined  by  the  intersection  of  lines 
in  the  xQyo  - plane  of  the  form  (21).  Of  course,  this  scheme  is 
prone  to  serious  numerical  errors  if  the  eddy  is  not  nearly 
circular,  if  the  constructed  lines  are  nearly  parallel,  or  if 
significant  errors  are  present  in  the  current  measurements. 

Effects  of  the  last  two  errors  can  be  avoided  by  using  several 
widely-spaced  current  meters,  and  avoiding  readings  where  |v| 
is  small.  In  practice,  measurements  of  currents  would  be 
discarded  if  the  current  speed  is  less  than  some  minimum  speed 
Vmin'  w^ere  vrain  ;''s  su^ f iciently  large  (perhaps  20  cm/sec)  that 
contributions  from  eddy  drift  or  other  environmental  effects 
are  insignificant. 

5.  Time- independent  problems 

Large  numbers  of  eddy  observations  consist  of  temperature 
measurements,  taken  over  sufficiently  short  time  intervals  so 
as  to  be  considered  time-independent.  Several  eddy  experiments 
have  made  use  of  ship-dropped  bathythermographs  (BTs)  to  get 
continuous  temperature  readings  to  depths  below  that  of  significant 
eddy  influence.  For  example,  much  data  from  the  MODE,  POLYGON, 
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and  POLYMODE  experiments  (see,  for  example,  Volkmann,  1977)  in 
the  form  of  temperature  sections  taken  in  linear  traversals  of 
an  ocean  area,  often  obtained  from  expendable  bathythermo- 
graphs (XBTs ) , is  available.  XBTs  are  less  expensive  and  easier 
to  use  than  BTs  but  are  restricted  to  depths  of  less  than  2 km. 
Also,  sections  are  obtained  from  airborne  expendable  bathythermo- 
graphs (AXBTs) , which  are  convenient  and  require  lower  support 
costs.  Although  they  are  restricted  to  depths  of  under  400  m 
at  present,  the  development  of  deeper  AXBTs  is  currently  being 
considered. 

We  divide  this  section  into  two  analytically  similar,  but 
experimentally  distinct,  parts.  First,  we  consider  data  from 
a vertical  planar  cross  section  through  an  eddy,  as  might  be 
obtained  by  a ship  linearly  traversing  an  eddy  or  by  a plane 
dropping  AXBTs  along  a linear  trajectory.  Then,  data  taken 
throughout  the  eddy  in  a nonplanar  fashion  is  analyzed.  The 
procedure  of  Secs.  3 and  4 and  the  model  of  Sec.  2 are  used 
to  approximate  the  corresponding  environmental  effects  of 

such  eddies. 

In  both  the  planar  and  nonplanar  cases,  the  data  is  assumed 
to  be  taken  over  sufficiently  short  time  intervals  that  eddy 
drift  may  be  neglected.  Thus,  the  coordinates  of  the  eddy 
center,  xQ(t)  = XQ,  yQ(t)  = Y , are  constant.  These,  together 
with  the  eddy  radius  rQ  and  maximum  current  speed  Uq,  result  in 
four  horizontal  parameters  (n  =4) , with  eddy  depth  of  influence 
zQ  as  the  sole  vertical  parameter  ( n v = 1).  The  direction  of 
rotation  is  chosen  immediately  by  the  sign  of  the  perturbation 
temperature. 

r 
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a.  Planar  temperature  data 

Sampling  of  oceanic  temperature  may  be  obtained  from 
ship-dropped  BTs  or  XBTs  or  plane-dropped  AXBTs  along  straight- 
line  traversals  of  large  ocean  areas.  Occasionally,  large 
cold-  or  warm-temperature  anomalies  are  observed,  whose  structure 
suggests  the  presence  of  a mesoscale  eddy.  Discrete  drops 
along  an  eddy  chord  are  illustrated  by  the  crosses  in  Fig.  la. 

A portion  of  an  XBT  section  through  the  Atlantic,  taken  by 
Seaver  (1975)  as  part  of  the  MODE  program,  is  shown  in  Fig.  2, 
indicating  the  presence  of  a large  cyclonic  eddy. 

As  proposed  in  Sec.  3,  zq  is  selected  first,  using  temperature 
measurements  at  the  same  horizontal  position  (x^,y^)  at  m^  (>2) 
depths  z^-^  , 1 £ j £ my.  Perturbation  density  is 

calculated  at  each  depth  as  in  (20) , and  is  used  as  the 
environmental  deviation  in  (15)  , with  the  theoretical  perturbation 
density  obtained  from  (18b) . The  solution  of  the  resulting  problem 

mv  f , 

Min  ^ |[AP;[j)-M  F*  (z{j))]2/(Ap{j))2}  (22) 

Zq,M  j=l 

then  specifies  the  vertical  eddy  structure  through  (7c)  and  (18) , 

At  each  of  the  m„  (>4)  horizontal  positions  (x.,y.),  and  depths 

ti  11 

zP^  , perturbation  density  is  calculated  as  in  (20).  Then  the 
average  vertically-scaled  perturbation  density  is  computed  from 
(13)  and  (18)  as 

"V 

A^i  = J [Ap|j)/F*  (z(j))]/mv  . 
j=l 


(23) 
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Determination  of  the  remaining  horizontal  parameters  follows  from 
(16),  using  Ap.^  from  (23)  as  the  average  vertically-scaled 
deviation,  resulting  in 


U ,r  ,X  , Y i=l 

o o o o 


£ { Api  -(2pong  sin  <t,)UoroKl[Jo(otlri/ro)"Jo(al)J 

i=1  — 1 

/(Ap.)2  J-  , (24a) 


where 


r.  = [ (x.-X  )2  + (y.-Y  )2J1/2 
1 1 o o 


(24b) 


Unique  determination  of  the  eddy  center  is  not  possible 

here.  As  indicated  in  Fig.  1(a),  two  circular  eddies  can  exist 

which  possess  the  same  temperature  structure  in  a common  vertical 

plane.  One  eddy  is  represented  by  the  solid  circle,  whose 

center  (X „,Y  ) is  a distance  d from  the  chord.  The  second, 
o o 

indicated  by  the  dashed  circle,  is  of  identical  size  and  strength, 

but  is  centered  at  (Xo,Yq)  at  a distance  d on  the  opposite  side 

of  the  chord.  Although  differentiation  between  the  two  centers 

is  not  possible,  values  for  the  other  horizontal  parameters, 

including  eddy  size  rQ  and  strength  Uq,  can  be  obtained.  Any 

convergent  scheme  used  in  solving  the  minimization  problems  will 

estimate  one  of  the  two  centers,  and  the  position  of  the  alternate 

center  can  be  found  by  reflection  across  the  eddy  chord.  The 

results  of  the  planar  temperature  data  problem  are  summarized 

in  the  first  line  of  Table  1.  In  this  example,  eddy  depth  of 

influence  z , radius  r , and  maximum  current  speed  U can  be 
o o o 
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uniquely  determined.  An  ambiguity  exists  in  the  time-independent 
position  of  the  eddy  center,  and  no  information  is  available 
concerning  the  drift  of  the  eddy. 

To  illustrate  our  procedure,  we  use  data  from  Fig.  2, 
taken  every  20  km.  At  each  horizontal  location,  the  temperature 
was  noted  at  seven  equally-spaced  depths  from  150  m to  750  m. 

Solving  (22)  numerically  gave  a depth  of  influence  of  zQ  = 2640  m. 
Eddy  size  and  strength  were  determined  by  solving  (24)  numerically , 
resulting  in  rQ  = 150  km  and  Uq  = 110  cm/sec.  The  distance  of 
closest  approach  of  the  chord  to  the  eddy  center  was  found  to  be 
10  km,  so  that  the  ambiguity  in  the  center  location  is  small  in 
this  example.  For  an  acceptable  horizontal  fit,  at  least  eight 
readings  were  necessary  to  obtain  physically  reasonable  parameter 
values.  Addition  of  more  points  led  to  better  agreement  in  the 
shape  and  elevation  of  the  isotherms,  although  inclusion  of  more 
than  twenty  points  led  to  insignificant  changes  in  parameter 
values.  An  example  of  the  dependence  of  parameter  values  on  the 
amount  of  data  will  be  presented  in  the  following  subsection. 
Resulting  isotherms,  obtained  from  (18),  the  ocean  static  state, 
and  the  inverted  Eckart  equation,  are  shown  in  Fig.  3,  Satisfactory 
agreement  is  noted  between  depths  and  overall  shapes  of  the  observed 
and  predicted  isotherms  in  Figs.  2 and  3.  Distributions  of  other 
eddy  characteristics,  such  as  currents,  can  be  determined 
similarly. 

b.  Non-planar  temperature  data 

We  consider  now  the  use  here  of  certain  non-planar 
temperature  data  from  BTs , XBTs,  or  AXBTs . Such  data  is  more 
useful  than  planar  measurements  in  experiments  designed  to  study 
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a particular  oceanic  eddy.  Many  previous  experimental  studies 
have  required  deep  temperature  data  to  determine  the  full 
vertical  structure  and  depths  of  significant  eddy  influence. 
However,  if  high  accuracy  at  large  depths  is  not  required,  then 
shallower  data  can  be  used  to  obtain  the  model  approximation  to 
the  vertical  eddy  structure.  With  the  model  of  Sec.  2,  we  shall 
focus  on  using  AXBTs,  although  the  results  will  be  valid  also 
for  BTs  or  XBTs. 

To  determine  the  vertical  structure,  we  have  noted  that 
measurements  to  depths  of  750  m are  necessary.  As  mentioned 
previously,  although  limited  to  depths  of  about  350  m at  present,  . 
the  possible  development  of  deeper  AXBTs  is  under  investigation. 
Thus,  estimation  of  zQ  by  solving  (22)  would  require  the  use  of 
one  such  AXBT,  or  another  form  of  measurement  of  the  eddy  vertical 
structure. 

With  the  eddy  vertical  structure  specified,  AXBTs  are  dropped 
at  (>4)  horizontal  positions  (x^,y^) , and  readings  are  obtained 
at  depths  z|^ . We  illustrate  in  Fig.  lb  one  scheme  for  obtaining 
nonplanar  data.  The  scheme  was  chosen  for  its  simplicity,  with  no 
attempt  being  made  here  to  select  an  optimum  procedure,  and  pre- 
supposes a nearly  circular  eddy.  A linear  path  L^,  with 
relatively  widely  spaced  AXBT  drops,  for  example  every  25  km,  is 
followed  until  several  temperature  anomalies  are  found,  indicating 
the  presence  of  a large  eddy.  Then,  readings  close  to  the  eddy 
center  are  obtained  by  making  more  closely  spaced  drops,  say 
every  10  km,  following  a path  L_  which  is  perpendicular  to  L. 

i 

and  passes  through  the  position  of  maximum  perturbed  temperature 

II 


■ 


28 


obtained  previously  along  L^.  The  second  chord  then  approximates 

a diameter.  Information  from  both  paths  permits  an  accurate 

determination  of  U , r , X , and  Y , by  first  computing  the 

o o o o 

average  vertically-scaled  perturbation  density  as  in  (20)  and 
(23)  at  each  drop  and  than  solving  problem  (24) . Moreover, 
ambiguity  in  eddy  center  position,  found  in  the  previous  sub- 
section, no  longer  exists.  Solution  of  the  horizontal  problem 
does  not  require  data  below  350  m,  and  hence  AXBTs  are  ideal  in 
providing  data  throughout  the  eddy,  after  the  vertical  structure 
is  specified.  Thus,  in  the  case  of  nonplanar  temperature  data, 
unique  determination  of  the  time-independent  parameters  is 
possible,  using  one  deep  measurement  at  a fixed  horizontal  position. 
We  summarize  these  results  in  the  second  line  of  Table  1,  Here, 
eddy  depth  of  influence,  radius,  maximum  current  speed,  and  position 
of  the  eddy  center  can  be  determined  at  the  time  of  the  measurements. 
However,  nothing  can  be  infered  concerning  the  translation  direction 
or  speed  of  the  eddy. 

As  an  illustrative  example,  we  use  BT  data  from  a large 
cyclonic  eddy  studied  by  Khedouri  and  Gemill  (1974) . The  eddy 
was  observed  to  have  significant  environmental  effects  to  depths 
of  over  2 km,  with  a radius  of  approximately  150  km,  and  a maximum 
current  of  110  cm/sec.  The  observed  temperature  data  in  a 
vertical  cross  section  through  the  eddy  center  was  rotated  to 
generate  an  axisymmetric  eddy.  A single  temperature  measurement 
(to  800  m)  near  the  eddy  center  was  used  to  determine  the  depth 
of  influence  The  resulting  value  of  2600  m is  in  agreement 

with  observations  of  the  vertical  extent  of  the  eddy.  The 
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equivalent  of  AXBT  data  can  be  constructed  by  sampling  the  data 
only  to  a depth  of  about  350  m.  Examples  of  typical  temperature 
deviations  at  several  distances  from  the  eddy  center  are  shown 
in  Fig.  4(a).  The  temperature  perturbations  increase  rapidly 
with  depth,  reaching  maxima  near  500  m (which  of  course  is  not 
observed  using  AXBTs) . A small  surface  expression  of  the  eddy, 
seen  in  the  observational  data,  and  presumably  due  to  near-surface 
mixing,  is  not  included  in  our  model  of  Sec.  2. 

Using  the  scheme  of  Fig.  1(b),  twenty-three  AXBTs  were 
simulated,  with  readings  taken  at  depths  of  150  m,  250  m,  and 
350  m,  at  the  positions  illustrated  in  Fig.  5(a).  We  obtained 
from  the  solution  of  (24)  the  values  UQ  = 120  cm/sec,  and  rQ  = 150  km. 
In  addition,  the  eddy  center  was  located  within  6 km  of  the  observed 
position.  Thus,  good  agreement  between  theory  and  observations 
was  obtained.  Theoretical  temperature  perturbations  are  shown  in 
Fig.  4(b).  They  exhibit  features  and  magnitudes  quite  similar  to 
those  observed  in  the  shallow  depths  of  Fig.  4 (a) , with  decay  at 
larger  depths  similar  to  that  observed  in  the  BT  readings.  In 
addition,  the  model  specifies  all  eddy  characteristics  at  all 
locations.  For  example,  Fig.  6 shows  curves  of  constant  rotational 
current  speed,  obtained  from  (18c) , in  a plane  through  the  eddy 
center.  The  currents  are  small  at  all  depths  near  the  eddy  edge 
and  center,  and  current  speeds  decay  rapidly  with  increasing  depth 
from  the  surface  maximum  at  all  horizontal  positions. 

The  number  of  points  necessary  for  accurate  parameter 
specification  was  tested  by  solving  the  horizontal  problem  (24) 
using  data  sets  of  increasing  size  taken  from  the  eddy  of  Fig.  4(a). 


30 


Initially,  three  AXBT  drops  were  made  on  LI  (Fig.  5(a))  and  two 

on  L2,  and  the  resulting  horizontal  problem  solved.  Additional 

AXBT  drops  were  made,  until  all  twenty-three  AXBTs  were  included. 

After  each  drop,  UQ,  rQ,  xq  and  yQ  were  determined  by  solving 

(24).  Resulting  values  of  UQ  and  rQ  are  plotted  in  Fig.  5(b) 

as  a function  of  the  number  of  AXBTs  utilized  in  the  horizontal 

fit.  In  this  example,  utilization  of  as  few  as  eight  AXBTs 

led  to  good  estimates  of  parameter  values.  However,  variations 

of  as  much  as  10%  in  parameter  values  from  the  "final"  values  of 

U = 120  cm/sec  and  r = 150  km  are  noted  if  fewer  than  twenty 
o o 

AXBTs  are  utilized.  Similar  results  are  observed  in  the  location 
of  the  eddy  center.  Of  course,  the  accuracy  of  the  parameter 
specifications  varies  not  only  with  the  number  of  AXBTs  used, 
but  also  with  the  positions  of  the  AXBTs,  the  eddy  under  con- 
sideration, and  the  accuracy  of  the  model  approximation.  We 
have  conducted  similar  studies  on  other  observed  eddies  and  with 
other  horizontal  locations  of  the  readings.  From  these  studies, 
we  have  found  that  a minimum  of  eight  AXBTs,  dropped  near  both 
the  eddy  edge  and  center,  are  ordinarily  necessary  for  good 
estimates  of  parameter  values.  Variation  of  parameter  values 
decreased  with  increasing  density  of  AXBTs,  until  inclusion  of 
more  than  twenty  AXBTs  led  to  insignificant  parameter  variation, 
as  in  Fig.  5 (b) . 

6.  Time-dependent  problems 

We  consider  here  the  analysis  of  eddy  observations  from 
fixed  moorings,  containing  instrument  packages  at  discrete  fixed 
depths.  In  contrast  to  BTs,  the  time  series  of  data  from  moorings 
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can  provide  information  on  eddy  translational  rates  and  directions, 
as  well  as  on  eddy  size  and  strength.  Temperature  observations 
are  available  from  thermistor  strings,  while  current-temperature 
moorings  provide  both  temperature  and  current  data;  we  will 
discuss  each  of  these  separately  in  Secs.  6(a)  and  6(b), 
respectively,  using  the  model  of  Sec.  2. 

a.  Thermistor-string  data 

We  first  consider  the  use  of  data  from  a single  thermistor 
string,  with  thermistors  located  at  (>2)  depths.  At  a fixed 
time  t^,  eddy  vertical  structure  may  be  determined  first  by 
solving  problem  (22).  At  each  time  t^,  the  average  vertically- 
scaled  perturbation  density  Ap^  is  then  determined  by  (23).  To 
further  describe  the  eddy,  we  would  like  to  specify  the  translation 
of  the  eddy  center  xQ(t)  and  yQ(t).  However,  using  a single 
thermistor  string,  it  will  be  possible  to  determine  only  the 
speed  by  which  the  eddy  drifts,  not  the  direction  of  translation. 
That  is,  only  the  eddy  drift  speed, 

2 2 X/2 

S(t)  = J (t)  | = {[x'(t)r  + [y ' (t)  ] Z } , (25) 

o o 

can  be  found.  We  shall  assume  here  that  the  drift  speed  is 
constant,  as  in  (17) , so  that 

2 2 

S = (u^  + Vp)  (26) 

is  one  horizontal  parameter  in  our  problem.  This  assumption  may 
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be  valid  if  the  eddy  is  observed  only  over  a sufficiently  short 
time  interval,  or  if  the  eddy  is  not  significantly  influenced  by 
other  ocean  currents  and  by  topographical  features.  Moreover, 
any  eddy  of  a fixed  strength  and  at  the  same  radial  distance  r 
from  the  thermistor  string  will  produce  identical  temperature 
effects.  Hence,  it  is  not  possible  to  determine  the  exact 
initial  position  (Xq,  Yq)  but  only  the  initial  distance  R from 
the  string.  Thus 

R = (X  2 + Y 2)1/2  (27) 

o o 

is  a second  horizontal  parameter.  Further,  only  the  distance  d 
of  closest  approach  of  the  eddy  center  to  the  string  can  be 
determined,  giving  a third  horizontal  parameter.  The  eddy 
radius  rQ  and  maximum  current  speed  Uq  comprise  the  remaining 
horizontal  parameters,  so  that  nH  = 5. 

To  determine  the  five  horizontal  parameters,  the  vertically- 

scaled  perturbation  density  Ap.  is  computed  at  m (>5)  times 

t. , i = l,...,m  . To  determine  the  predicted  effects  through 
X n 

(18b) , the  radial  distance  r must  be  computed  at  each  time.  We 
define  a moving  coordinate  system  (x',y'),  with  origin  at  the  eddy 
center.  The  x'  -axis  is  parallel  to  the  direction  of  eddy  drift, 
with  decreasing  x'  in  the  direction  of  Vp,  as  illustrated  in 
Fig.  7(a).  To  an  observer  fixed  in  the  (x',y')  system,  the 
thermistor  string  will  appear  to  translate  in  the  positive  x' -direction 
with  speed  S.  Initially,  the  eddy  center  will  be  located  at 
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(Xo',Yo')  = [-(R2-d2)1/2,-d]  . (28a) 

From  elementary  geometry,  the  desired  radial  distance  r(t)  at 
any  time  t will  then  be 

1/2 

r (t)  = [(Xo'+St)2  + d2]  . (28b) 

Equations  (28)  together  with  (24) , comprise  the  problem  for  the 
five  horizontal  parameters.  Examination  of  the  problem  shows 
that  it  can  be  solved  only  for  the  parameter  combinations 
UQro,R/ro,d/ro,  and  S/rQ.  That  is,  unique  solutions  for 
individual  parameter  values  are  not  possible,  since  UQ,R,d, 
and  S can  be  determined  only  as  (simple)  functions  of  tq.  Thus,  . 
each  member  of  a one-parameter  family  of  eddies  is  a possible 
fit  to  the  observed  temperature-perturbation  data.  The  nonunique 
parameter  specification  is  summarized  in  the  third  line  of  Table  1. 
Eddy  depth  of  influence  can  be  determined  uniquely,  but  eddy 
radius,  maximum  current  speed,  and  the  eddy  center  position  and 
trajectory  can  only  be  determined  to  within  a one-parameter 
family  of  values,  assuming  a linear  drift  form. 

To  illustrate  this  nonuniqueness,  we  present  simulated 
thermistor  string  data  in  Fig.  8 for  a large  cyclonic  eddy  with 
linear  drift.  The  data  for  the  resulting  cold-core  eddy  was 
generated  using  the  model  of  Sec.  2,  assuming  a drift  speed 
S = 5 km/day,  radius  rQ  = 125  km,  maximum  current  speed  UQ  = 150 
cm/sec,  closest  approach  distance  d = 20  km,  and  depth  of  influence 
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zq  = 2100  m.  We  note  that,  at  any  depth,  temperature  decreases 

as  the  eddy  approaches  the  thermistor  string,  and  then  increases 

to  the  static  state  as  the  eddy  recedes.  At  the  initial  time 

t = 0,  the  thermistor  was  chosen  to  be  on  the  eddy  edge,  so  that 

R = rQ.  Applying  the  procedure  of  this  section,  readings  taken 

at  the  depths  shown  on  Fig.  8 every  five  days  reproduce  the 

exact  (inputed)  value  of  zQ  = 2100  m,  and  yield  constant  values 

for  U r ,d/r  , and  S/r  . The  parameters  U ,d,  and  S are 
o o o o o 

illustrated  as  functions  of  rQ  in  Fig.  9.  We  note  that  for  any 
assigned  '^’ue  of  the  eddy  radius,  S,  d,  and  Uq  are  uniquely 
specified,  as  indicated  by  the  dashed  line  corresponding  to  the 
exact  radius  rQ  = 125  km.  To  generate  curves  such  as  those  in 
Fig.  9,  a minimum  of  four  readings  at  different  times  is  required; 
for  accurate  specification,  we  have  found  that  at  least  eight 
readings  are  necessary,  each  separated  by  several  days. 

Addition  of  more  than  one  thermistor  string  at  different 
horizontal  positions  supplies  sufficient  information  to  determine 
unique  parameter  values,  and  to  specify  drift  direction  as  well 
as  speed.  The  vertical  parameter  zq  is  again  chosen  first  by 
using  data  from  any  one  string  at  a fixed  time.  I*  we  assume  the 
linear  trajectory  (17),  the  radial  distance  r^ (t)  from  the 
thermistor  string  T^  to  the  eddy  center  at  time  t is  clearly 

f 2 2 1 

ri(t)  = |[xi  - (Xo  + uDt)p  + [y.  - <Yo  + vDt)  ] J . (29) 

To  determine  the  horizontal  parameters  UQ,  rQ,  Xq,  Yq,  ud,  and 


V n»H  (>6)  readings  are  required.  Using  measurements  from  at 
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least  two  thermistor  strings,  a unique  solution  can  be  obtained, 
since  parameters  now  occur  in  (24a)  and  (29)  in  the  six  distinct 
combinations  UQro,  l/rQ,  Xo/ro,  Y q/tq.  uD/rQ,  and  vD/rQ.  To 
accurately  specify  model  parameters,  a significant  overlap  in 
time  (say  15  days)  from  data  taken  from  different  strings  must 
be  present,  and  separation  between  thermistor  strings  should 
not  be  small  in  comparison  to  the  size  of  the  eddy  (for  example, 
50  km).  Thus,  by  the  inclusion  of  additional  thermistor  strings, 
the  indeterminacy  of  the  one  string  case  is  avoided,  and  unique 
parameter  specification  is  possible.  Results  are  summarized 
in  the  fourth  line  of  Table  1.  Eddy  depth  of  influence,  radius, 
maximum  current  speed,  center  position  and  drift,  assuming  linear 
drift,  are  uniquely  determined.  More  complex  forms  for  the  eddy 
translation  xQ(t)  and  yQ (t)  may  be  assumed,  leading  to  additional 
parameters,  which  can  also  be  uniquely  determined  using  two  or 
more  thermistor  strings. 


Temperature  and  current  data 


Current  data  provides  significantly  more  information  than 
is  available  from  temperature  measurements  alone.  The  inclusion 
of  current  meters  on  a thermistor  string  provides  additional 
information  on  both  the  eddy  strength  at  a given  time,  through 
current  speed  observations,  and  the  location  of  the  eddy  center, 
via  current  direction. 

We  first  consider  a single  temperature-current  mooring  with 

n»c  current  meters  at  depths  n ^ 1 £ 3 £ mc>  anc*  ^ thermistors 
( o ) 

at  depths  z , 1 £ Jl  <_  n^.  At  time  t^,  the  current  direction 
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and  magnitude  at  the  j*"*1  meter  are  and  The  vertical 

structure  can  again  be  chosen  first  by  using  temperature 
readings  at  a fixed  time  t^  and  solving  problem  (22).  We  note 
that  current  magnitude  could  also  be  used  as  the  environmental 
deviation  in  such  a determination,  in  which  case  the  problem 
becomes 


Min 


[v|j)  - MlF(n(j))-F(zo)]j2/(V^j))2 

j=l 


(30) 


where  F(z)  is  given  by  (7c).  We  caution  again  against  using 
readings  where  current  magnitudes  are  less  than  Vm^n  (say  20  cm/sec) , 
as  errors  induced  by  eddy  drift  and  other  environmental  effects 
might  become  significant.  Alternatively,  both  current  and 
temperature  data  could  be  used  in  selecting  zq,  minimizing  the 
sum  of  the  squares  of  both  (22)  and  (30) . At  any  time,  the 
average  vertically-scaled  perturbation  density  can  be  determined 
and,  similarly,  the  average  vertically-scaled  current  speed 
given  as 

m 

c 

V . = (1/m  ) Y V.(j)/[F(n(j,)-F(z  )]  (31) 

1 v*  * 1 O 

j=l 

where  readings  with  currents  below  V . are  not  included. 

min 

Current  direction  measurements  can  then  be  used  to  simplify 
the  horizontal  problem.  At  any  time  t^,  the  average  direction 
6^  can  be  obtained  by  averaging  the  current  direction  from  each 
meter  for  which  . The  relationship  (21)  gives  a line 
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of  possible  centers,  giving  yQ(t  ) in  terms  of  x^t^,  and 

hence  eliminating  the  parameters  implicit  in  yQ(t).  For  the 

form  of  x (t) , we  again  assume  a linear  drift,  (17),  although 
o 

a nonlinear  form  can  be  treated  by  introducing  the  appropriate 

parameters.  Thus,  there  are  four  horizontal  parameters, 

X ,u  ,r  , and  U , and  the  sole  vertical  parameter  z . From 
o D o o o 

(21),  the  radial  distance  from  the  eddy  center  at  time  t^  is 
given  by 

r(ti)  = |XG  + uDti|  [1  + tan2  (ir/2  - Q±)}1/2.  (32) 


We  note  from  Fig.  6 that,  near  the  eddy  center  or  edges,  eddy 
currents  at  all  depths  are  small.  Hence  the  use  of  current 
direction  to  eliminate  the  parameters  in  yQ(t),  through  (32), 
should  be  avoided  in  these  regions.  The  problem  for  determining 
the  remaining  model  parameters,  given  readings  at  mH  (>4) 
distinct  times  when  significant  currents  (>V  . ) are  present  is 


. ,u  ,r  ,U  Zl 
o D O O 


f J(A0i  - (20ofi  9_lsin  4>)  Uoro  VVV7 r0) 


-Jo(ol)1)2/(lpi|2  + I?i  ‘ "o  Kl“lJl(t>lr''ro)]2/'?i} 


where  r(t^)  is  obtained  from  (32).  Since  UQ  now  occurs  distinctly 
in  the  current  term,  unique  solutions  can  be  found.  This  was  not 
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possible  when  only  temperature  data  was  available.  Eddy  depth 
of  influence,  radius,  maximum  current  speed,  center  position 
and  linear  drift  can  all  be  determined  uniquely,  as  summarized 
in  line  5 of  Table  1. 

Suppose  now  that  n (>1)  current-temperature  moorings  are 
present  at  any  time  t.  From  the  discussion  of  Sec.  4,  the 

position  of  the  eddy  center  (xq (t) ,yQ  (t) ) can  be  determined  by  1 

the  intersection  of  perpendiculars  to  the  average  vertically- 
scaled  current  at  each  mooring  computed  as  in  (31) . This  is 
indicated  in  Fig.  7(b)  for  the  case  of  two  current  meters. 

The  position  of  the  eddy  center  is  known  whenever  two  or  more 
current  meters  measure  significant  currents;  for  accurate 
specification  of  the  center  position,  readings  should  be  taken 
from  more  than  two  moorings  at  several  times.  Of  course, 
because  of  observational  inaccuracies,  if  n > 2,  the  (n-l)i 
intersection  points  would  not  coincide.  They  would  be  expected 
to  be  close  enough  so  that  a reasonable  estimate  of  the  center 
could  be  made,  if  the  guidelines  of  Sec.  4 are  followed; 
otherwise,  the  radially  symettric  model  may  not  be  applicable 
to  the  observed  eddy.  The  vertical  structure  is  specified  next 
by  using  readings  from  a single  mooring  at  a fixed  time,  and 
solving  problem  (22)  or  (30) . Average  vertically-scaled 
perturbations  and  current  speeds  are  then  computed  at  each 
mooring  and  time,  and  a problem  similar  to  (33)  solved  for  the 
remaining  horizontal  parameters  UQ  and  rQ.  Thus,  with  two  or 
more  current  moorings,  the  eddy  position  at  any  time,  as  well 
as  the  size  and  strength,  can  be  determined.  This  is  summarized 
in  the  last  line  of  Table  1.  Eddy  depth  of  influence,  radius. 
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maximum  surface  current  speed,  and  eddy  center  position  and 
drift  can  all  be  determined,  independent  of  any  assumptions 
concerning  the  form  of  the  drift  trajectory. 

7 . Summary 

The  major  purpose  of  this  paper  is  to  describe  the  use 
of  analytic  models  to  determine  approximate  environmental 
effects  of  an  oceanic  eddy  using  limited  observational  data. 

The  question  of  the  amount  and  types  of  data  necessary  for  such 
model  specification  is  discussed,  and  techniques  for  accurate 
model  parameter  specification  presented. 

In  order  to  provide  illustrative  examples  of  the  ideas 
in  this  paper,  a previously-derived  model  that  gives  predictions 
of  eddy  currents,  density,  pressure,  and  temperature  is  briefly 
reviewed.  This  model  has  eddy  radius,  maximum  current  speed, 
depth,  and  drift  trajectory  as  parameters , but  any  model  with 
an  arbitrary  number  of  parameters  can  be  utilized.  Thus, 
parameter  specification  for  a general  model  is  discussed. 
Separable  models  are  considered,  in  which  parameters  associated 
with  the  vertical  and  horizontal  structure  are  determined 
separately.  The  minimum  amount  of  data  required  to  specify 
an  eddy  model,  a necessary  prerequisite  to  accurate  model 
approximation,  is  analyzed.  Then,  an  efficient  strategy  is 
developed  for  accurate  determination  of  model  parameters, 
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employing  successive  minimization  problems  for  the  vertical 
and  horizontal  parameters.  These  parameters  and  the  equations 
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specifying  them  are  then  described  for  the  eddy  model  previously 
presented. 

The  use  of  oceanographic  data  in  parameter-specification 
schemes  is  discussed.  Although  the  direct  use  of  temperature 
data  can  result  in  intolerable  numerical  errors,  density 
anomalies  induced  by  an  eddy  can  be  used  to  accurately  determine 
model  parameters.  The  use  of  current  measurements  are  discussed 
also.  The  procedures  are  then  applied  to  several  typical 
experimental  situations,  using  our  particular  eddy  model  as  an 
example.  First,  temperature  data  acquired  on  a linear  path 
through  the  eddy,  as  might  be  obtained  during  a single  traversal 
of  an  ocean  region,  is  considered.  Assuming  that  the  data  are 
taken  over  small  enough  intervals  to  be  time-independent,  it  is 
shown  that  the  size  and  strength  of  a particular  eddy  can  be 
determined  uniquely.  However,  an  ambiguity  results  in  the 
position  of  the  eddy  center,  so  that  either  of  two  possible  eddies 
could  be  responsible  for  the  observed  temperature  perturbations. 
Next,  we  study  non-planar  time-independent  data,  with  emphasis 
on  that  from  AXBTs , as  might  be  obtained  during  an  extensive 
study  of  an  eddy.  It  is  shown  that  unique  determination  of  eddy 
size,  strength,  and  position  is  possible,  although  no  information 
may  be  learned  about  the  eddy  drift.  Examples  are  presented  from 
actual  eddy  observations,  with  eddy  size  and  strength  determined 
to  within  five  percent  of  observed  values,  and  center  position 
accurate  to  within  10  km.  Guidelines  are  given  for  placement  of 
bathythermographs  for  accurate  parameter  specification,  and 
numerical  sensitivity  to  the  amount  and  position  of  data  is 
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considered. 

Two  distinct  types  of  time-series  data  are  examined,  the 
first  being  only  temperature  readings  that  might  be  obtained 
from  one  or  more  thermistor  strings,  and  the  second  including 
both  temperature  and  current  readings  from  one  or  more  moorings. 

For  a single  thermistor  string,  it  is  shown  that  the  model 
parameters  cannot  be  determined  uniquely,  but  that  possibilities 
are  restricted  to  a one-parameter  family  of  eddies.  This 
includes  drift  speed,  but  not  direction,  if  a constant  drift 
velocity  is  assumed.  Addition  of  more  strings,  however,  does 
lead  to  unique  parameter  values,  as  well  as  specification  of 
drift  direction.  Addition  of  current  meters  to  a single 
thermistor  string  leads  to  complete  eddy  specification  when 
linear  drift  is  assumed.  Finally,  data  from  two  or  more  such 
moorings  can  be  used  to  determine  arbitrary  drift  as  well  as 
eddy  depth,  radius,  and  maximum  rotational  current  speed. 

The  use  of  analytical  eddy  modeling  simplifies  the  amounts 
and  types  of  data  required  to  determine  eddy  size,  strength, 
and  motion,  and  permits  efficient  approximating  of  eddy  environ- 
mental effects.  Future  work  in  refining  the  parameter  specification 
technique  should  be  in  the  areas  of  model  improvement  and  testing 
with  other  eddy  models,  further  sensitivity  analysis  by  comparison 
of  predictions  and  observations  using  different  and  varied  data 
sets  and  estimation  procedures,  and  utilization  of  the  technique 
with  other  forms  of  observations  such  as  float  paths.  Hopefully, 
such  refinements  will  contribute  to  the  use  of  analytical  modeling 
in  describing  ocean  processes. 
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LEGENDS 

(a)  Chord  through  eddy,  showing  ambiguity  in  eddy 
center  locations,  (b)  one  scheme  for  obteiining  non- 
planar  eddy  temperature  data. 

Experimental  isotherms  (interpolated)  through  a 
large  cyclonic  eddy  (adapted  from  Seaver,  1975). 

Theoretical  isotherms  for  an  eddy  with  radius  150  km, 
maximum  current  speed  110  cm/sec,  and  depth  of 
influence  2640  m. 

Eddy  temperature  perturbations:  (a)  "observational" 
AXBT  data,  (b)  theoretical  results. 

(a)  Position  of  AXBT  drops  in  eddy  of  Fig.  4,  (b) 

variation  of  eddy  radius  and  maximum  current  speed 
with  increasing  numbers  of  AXBTs. 

Theoretical  current  structure  through  an  eddy  center 
with  radius  147  km,  maximum  current  speed  123  cm/sec, 
and  depth  of  influence  2600  m. 

(a)  One  thermistor  string  geometry,  (b)  location  of 
eddy  center  using  two  current  meters. 

Model-simulated  thermistor  string  data  from  eddy  with 
linear  drift  speed  5 km/day,  radius  125  km,  maximum 
current  speed  150  cm/sec,  closest-approach  distance 
25  km,  and  depth  of  influence  2100  m. 

Family  of  parameter  values  for  eddy  of  Fig.  7. 
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